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ABSTRACT 

We numerically study the volume density probability distribution function 
(ra-PDF) and the column density probability distribution function (S-PDF) re- 
sulting from thermally bistable turbulent flows. We analyze three-dimensional 
hydrodynamic models in periodic boxes of 100 pc by side, where turbulence is 
driven in the Fourier space at a wavenumber corresponding to 50 pc. At low den- 
sities {n < 0.6 cm~^) the n-PDF, is well described by a lognormal distribution 
for average local Mach number ranging from ~0.2 to ~5.5. As a consequence 
of the non linear development of thermal instability (TI), the logarithmic vari- 
ance of the distribution for the diffuse gas increases with M faster than in the 
well known isothermal case. The average local Mach number for the dense gas 
{n > 7.1 cm~^) goes from ~1.1 to ~16.9 and the shape of the high density zone 
of the n-PDF changes from a power-law at low Mach numbers to a lognormal at 
high M values. In the latter case the width of the distribution is smaller than 
in the isothermal case and grows slower with M. At high column densities the 
E-PDF is well described by a lognormal for all the Mach numbers we consider 
and, due to the presence of TI, the width of the distribution is systematically 
larger than in the isothermal case but follows a qualitatively similar behavior as 
M increases. Although a relationship between the width of the distribution and 
M can be found for each one of the cases mentioned above, these relations are 
different form those of the isothermal case. 
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1. 



Introduction 



The density distribution and particularly the density Probability Distribution Function 
(PDF) has become a crucial ingredient in theories about molecular cloud and core formation 
and evolution and on star formation theories (e.g. Padoan Sz Nordlund 2002; Krumholz & 
McKee 2005; Elmegreen 2002, 2008; Hennebelle & Chabrier 2008; Elmegreen 2011; Padoan 
& Nordlund 2011; Zamora et al. 2012). 

The volume density PDF in turbulent compressible flows has been widely studied for 
the isothermal case. In particular, numerical experiments of driven and decaying turbulence 

(Vazquez-Semadeni 1994; Padoan ct al. 1997; Stone et al. 1998; Klessen 2000; Ostriker 
et al.2001; Boldryev ct al. 2002; Bcetz et al. 2008; Federrrath ct al. 2008) have shown 
that this distribution has a lognormal shape in a large number of situations. From the 
theoretical point of view, the development of a lognormal distribution for isothermal flows 
have been explained as a consequence of the multiplicative central limit theorem assum- 
ing that individual density perturbations are independent and random (Vazquez-Semadeni 
1994; Passot & Vazquez-Semadeni 1998; Nordlund & Padoan 1999). Furthermore, Passot & 
Vazquez-Semadeni (1998) showed that the density PDF develops a power-law tail at high 
(low) densities for values of the polytropic index 7, smaller (larger) than 1. Also for the 
isothermal case, numerical simulations at a fixed rms Mach number M, show an empirical 
relationship between the width of the density PDF, as measured by its variance or its stan- 
dard deviation a and M (e.g. Padoan et al. 1997, Passot & Vazquez-Semadeni 1998). For 
the distribution of Inn/n this relationship reads 



where 6 is a constant of the order of unity whose value is not clearly established. In the 
literature it goes from 0.26 (Kritsuk et al. 2007) to 1 (Passot Sz Vazquez-Semadeni 1998) 
and seems to vary depending on the relative degree of compressible and solenoidal modes 
of the forcing (Federrath et al. 2008, Federrath et al. 2010). Recently Price et al. (2011) 
found that for M up to 20 6 = 1/3 is a good fit to numerical simulations with solenoidal 
forcing and is in agreement with other recent numerical results (Federrath et al. 2008, 2010; 
Lemaster & Stone 2008). 

On the other hand, the condensation of diffuse gas produced by the isobaric mode of 
thermal instability (TI; Field 1965) triggered in coUiding flows has been recognized as playing 




ln(l + 6'M2), 



(1) 
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an important role for molecular cloud formation (e.g. Hennebelle & Perault 1999; Vazquez- 
Scmadcni ct al. 2006; Heitsch & Hartman 2007), specially in early phases. However the 
density PDF resulting in turbulent bi-stable flows, which preserves the bi-modal nature that 
is signature of the development of TI, and the behavior of the distribution for each phase 
has not been systematically studied. 

The density PDF resulting from non-isothermal simulations, has been nevertheless re- 
ported in some papers. Using two-dimensional simulations of turbulent thermally bistable 
flows Gazol et al.(2005) have shown that the effective polytropic index of the gas increases 
with M, and that for Mach numbers between 0.5 and 1.25 (with respect to the gas at 10*^ K) 
it remains < 1, with specific values depending also on the forcing scale going from f» 0.2 to 
0.4 for large scale forcing. They also found that the bimodal nature of the distribution 
becomes less pronounced as the M increases. The same behavior has been reported from 
three dimensional simulations (Seifried et al. 2011). Audit & Hennebelle (2010) compare 
the density distribution resulting from a bi-stable simulation and the one resulting when a 
single polytropic equation of state with 7 = 0.7 is used. In the latter case they find that 
the low density part of the PDF is well described by a lognormal distribution while for the 
higher densities the PDF is a power-law with an exponent of about -1.5, which is consistent 
with the distribution expected for a very compressible gas with 7 < 1 predicted by Passot 
& Vazquez-Semadeni (1998). For the coohng run they find that the density distribution of 
the cold gas, whose effective polytropic index is close to 0.7, is well fitted by a log-normal 
for values larger than ~ 300 cm~^. They attribute the difference in the behavior to possi- 
ble resolution effects and rise the question about whether or not the choice of log-normal 
distribution in models of molecular clouds is adequate. 

The density distribution resulting from more complex systems including a variety of 
physical ingredients has been studied by a large number of authors. Here we just mention 
some examples that illustrate the activity in this direction. In the context of galactic discs, 
the effect of varying the polytropic index has been numerically studied by Li et al. (2003) 
who find that for self-gravitating gas the density distribution shows an imperfect lognormal 
shape whose width decreases as 7 increases. Wada & Norman (2007) use three-dimensional 
hydrodynamic simulations of a globally stable, inhomogeneous cooling ISM in galactic disks 
and report density PDFs which are well fitted by a single lognormal function (with larger 
dispersions for more gas-rich systems) over a wide density range. However, Robertson & 
Kravtsov (2008), include a detailed model of cooling and heating for temperatures < 10^ K 
and account for the equilibrium abundance of H2. They conclude that each thermal phase 
in their model galaxies has its own lognormal density distribution, implying that using 
a single lognormal pdf to build a model of global star formation in galaxies is likely an 
oversimplification. Density distributions with well defined peaks have also been obtained 
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from models of vertically stratified interstellar medium dominated by supernovae (e.g. De 
Avilez & Breitschwerdt 2005; Joung & MacLow 2006). 

Observationally, PDFs of the average volume density in diffuse interstellar gas (nni < 
1 cm~^ for the atomic gas) have been reported by (Berkhuijsen & Fletcher 2008). For the 
atomic gas they used 375 lines of sight from the sample given by Diplas & Savage (1994) and 
they estimate the volume density using the column density and the distance to each star. 
They find that the PDF tends to be a log-normal but cannot be well fitted by a single curve. 
In particular different widths and means are obtained depending on the position (in the disk 
or away from it) and on the type of sampled gas. Specifically, they find that a mixture of 
cool and warm gas along the LOSs causes an increase in the dispersion. 

An associate tool is the column density PDF (S-PDF), which has been numerically 
studied mainly in the context of isothermal turbulent fiows. Vazquez-Semadeni & Garcia 
(2001) studied the relationship between a lognormal volume density PDF (n-PDF) and the 
resulting column density PDF, finding that when the number of decorrelated density struc- 
tures along the line of sight is small enough, then the S-PDF is representative of the n-PDF 
and do also have a lognormal shape. As the number of decorrelated density structures in- 
creases the column density PDF slowly transits to a Gaussian distribution passing through 
an intermediate stage where the distribution shows an exponential decay. As the relation- 
ship between the ?7,-PDF and the S-PDF depends on the applicability of the Central Limit 
Theorem, the authors argue that for non-isothermal fiows having a volume density PDF 
with a well defined variance, this theorem should apply and the S-PDF should converge to 
a Gaussian for lines of sight with a large enough number of decorrelated density structures. 
More recently, Burkhart & Lazarian (2012) used solenoidaly driven isothermal MHD simu- 
lations to investigate the presence of an empirical relationship between the variance of the 
column density distribution and the sonic Mach number. They found a relationship with 
the same form as equation ([3]), namely 

<E/Eo) = ^sMl + &|M^), (2) 

where the scaling parameter = 0.11 and = 1/3, which is also close to the value 
reported for the 3D density distribution in the isothermal case . 

On the other hand, Ballesteros-Paredes et al. (2011) numerically studied the evolution 
of the S-PDF resulting from colliding self-gravitating fiows in presence of TI, which are in- 
tended to model the formation process of molecular clouds. They found that a very narrow 
lognormal regime appears when the cloud is being assembled. However, as the global grav- 
itational contraction occurs, the initial density fiuctuations are enhanced, resulting, first, 
in a wider lognormal S-PDF, and later, in a power-law S-PDF. These results suggest an 
explanation to the observational fact that clouds without star formation seem to possess a 
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lognormal distribution, while clouds with active star formation develop a power-law tail at 
high column densities (e.g Kainulainen et al. 2009). 

In this work we quantitatively study the behavior of the density and the column density 
distributions resulting in thermally bistable turbulent flows. For this purpose we analyze 
simple numerical 3-dimensional experiments which include only a cooling function appro- 
priate for the diffuse neutral interstellar gas and turbulent forcing. The paper is organized 
as follows: The model we use is described in section §21 Then in section §3] we present the 
analysis of the behavior of both the volume density PDF and the column density PDF. These 
results are discussed in section ^ Finally in ^we present our conclusions. 

2. The model 

We use the same model as in Gazol &; Kim (2010), where a MUSCL-type scheme (Mono- 
tone Upstream-centered Scheme for Conservation Laws) with HLL Riemann solvers (Harten, 
Lax, & van Leer 1983; Toro 1999) is employed to solve hydrodynamic equations in three di- 
mensions within a cubic computational domain with a physical scale of 100 pc by side. The 
turbulence is randomly driven in Fourier space at large scales corresponding to 1 < /c < 2, 
where k is the magnitude of the wave vector, k. We use a purely solenoidal forcing because 
of two reasons. First, it is the most common kind of forcing used in the isothermal case 
and second, also for the isothermal case is the kind of forcing for which the density PDF 
is better described by a lognormal distribution (Federrath et al. 2008). For further details 
concerning the turbulent forcing see Gazol & Kim (2010) and for a discussion on the effects 
of forcing see §4.4.31 In the real ISM, turbulent motions are driven by spatially localized 
sources, as the forcing in Fourier space is applied at any time in all the grid points, this 
kind of driving is not very realistic. However we chose this kind of driving method because 
it allows us to compare our results with previous work done for the isothermal case which 
includes Fourier space forcing. As an additional advantage, with this method we can study 
the density distribution for the dense gas as well as the density distribution for the diffuse 
gas (see We utilize the radiative cooling function presented by Sanchez- Salcedo et al. 
(2002), which is based on the standard P vs. p curve of Wolfire et al. (1995). 

For all simulations we present in the next section the resolution is 512^, the boundary 
conditions are periodic, the gas is initially at the rest, and the initial density and temperature 
are uniform with no = 1 cm~^ and Tq = 2399 K, which correspond to thermally unstable gas. 
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3. Results 

In this section we analyze the density distribution resulting from one set of seven sim- 
ulations with different Mach number. The values of M for those simulations are 0.28, 0.73, 
1.86, 2.60, 3.29, 5.50, and 6.23, where M is computed as the mean value of the local Mach 
number. In what follows we call this value (M) and we use M as a generic abbreviation for 
Mach number in situations that are independent of the specific way used to compute it or 
in situations where Mach numbers calculated in different ways are included. 

3.1. The Volume Density Distribution 

Volume density histograms resulting from these simulations are displayed in Figure [T] 
along with lognormal fits to the high and the low density parts of the distribution. As 
expected from previous works (e.g. Gazol et al. 2005; Seifried et al. 2011), when the Mach 
number increases the PDF becomes wider and its bimodal nature, which is a consequence 
of Tl development, becomes less pronounced. These two facts have been explained as a 
consequence of the decrease in the local ratio between the turbulent crossing time and the 
cooling time rj produced by the increase of M (Gazol et al. 2005). This behavior have 
been proved by Seifried et al. (2011), who using test particles measured the time spent 
by a particle in the unstable regime as well as the frequency with which a test particle is 
perturbed, finding that both quantities increase with M (i.e. both quantities increase as rj 
decreases). Concerning the shape of each part of the distribution, there are four things to 
note by a simple inspection. First, the low density zone of the PDF can be relatively well 
described by a lognormal (see dotted lines in Figure [1]). Second, the width of this lognormal 
increases with (M). Third, for the smallest values of (M) we consider, the zone of the 
distribution corresponding to high densities is not well described by a lognormal (see dashed 
lines in Figured]). Fourth, when a lognormal is an appropriate fit to the high density part, 
its width does not seem to systematically increase with the value of M as rapidly as in the 
diffuse case. 

In Figure [2] we show the width of the lognormal fits ag (note that cXs is the width of 
the logarithmic density distribution Inn/n) that we obtain for the gas at low densities as a 
function of the local Mach M„, which is computed as the average value of the local Mach 
number in points with densities within the range where the fit has been computed. The 
dotted lines in this figure are plots of the relationship, between cTs and M, for isothermal gas 
(see eqn. [T]) for two different values of the parameter h. Although for low values, is 
close to the 6 = 1/2 isothermal case, for our simulations it grows faster with M^. In fact. 
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red line represents a fit with the form 

aj = Aln{l + b^M^), (3) 

where A = 2.25 and b = 0.33. This value of b is close to typical values found for isothermal 
turbulence with solenoidal forcing. 

For the dense gas the behavior of cr^ with (computed as the average value of the 
local Mach number in points with n > 7.1 cm~'^, this density corresponds to the thermal 
equilibrium value below which the cooling function implies that the gas is thermally unsta- 
ble), is qualitatively different from equation ([3]) (see Fig. [3]). In this case, the distribution is 
narrower than in the isothermal case and its width grows more slowly with than in the 
isothermal case. In fact, due to the large difference between the widths we obtain and the 
ones implied by equation (|3]) we do not include both cases in the figure. For instance, for 
M = 10 in the isothermal case a.^ = 1.58 is expected if 6 = 1/3. 

As mentioned in ^11 previous works have found that the shape of the density PDF 
in non-isothermal turbulent flows depends on the value of the effective polytropic index 7 
(Passot & Vazquez Semadeni, 1998). For our simulations we first measured this index as a 
line least squares slope of the whole log P vs. log n distribution and we plotted it as a function 
of (M) (Fig. m solid line). As expected from Gazol et al. (2005), 7 increases with (M) and 
remains < 1. As can be seen from the log P/k vs. logn distributions displayed in Figure 
for low values of (M) TI can develop almost without disturbance and the dense as well as the 
diffuse gas are predominately in thermal equilibrium with some of it transiting isobarically 
between stable branches. On the other hand, at high (M) the mean pressure in low density 
gas decreases below the thermal equilibrium value while the mean pressure in high density 
gas increases above its thermal equilibrium value, producing a neat positive slope of the 
whole log P/k vs. logn distribution. If the value of 7 obtained in this way were a suitable 
parameter to infer the shape of the density PDF, we would expect a power law at high density 
for all our simulations. Also from Figure [5]it is clear that for low values of (M) almost all the 
dense gas {n > 7.1 cm~^) is in thermal equilibrium, implying that its thermodynamic state 
is approximately well described by a polytropic relation with 7 = 0.53, which is the power 
corresponding to the thermal equilibrium curve and which accordingly with the Passot & 
Vazquez-Semadeni (1998) theory should exhibit a power-law PDF. The dotted line in Figure 
m corresponds to the slope of the log P/ k vs. log n distribution computed for densities larger 
than 7.1 cm~^ and plotted as a function of the average local Mach number for dense gas 
Mc- This slope shows a decrease with from ~ 0.53 for (M) = 0.28 (Mc = 1.07) to 
~ for (M) = 6.23 (Mc = 16.91). This decrease is due to the increased dispersion at 
the lower density values of high density gas, and implies that the polytropic description 
of the dense gas becomes less appropriate for large M values. For diffuse gas the slope 
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of the distribution {dashed line), estimated for each simulation using density values within 
the range for which the fit to the corresponding PDF has been computed, increases from 
its thermal equilibrium value, namely 0.73, for low simulations to approximately 1 for 
very supersonic simulations, implying that the approximately lognormal shape we obtain is 
consistent with the predictions done by Passot & Vazquez-Semadeni (1998). 



3.2. The Column Density Distribution 

As long as the authors know, the column density distribution resulting from turbulent 
thermally unstable simulations, has not been systematically studied. It is thus interesting 
to study the S-PDF associated with the n-PDFs discussed in previous section. 

In Figure |6] we show the column density distribution resulting from our set of simulations. 
The high density part of the distribution can be well described by a lognormal {dotted lines) 
for all the values of (M) that we consider, even for the two smaller values (0.28 and 0.73) 
for which de volume density PDF has a behavior approaching a power-law. For these two 
values the bimodal nature of the n-PDF is preserved in the S-PDF. For larger values of 
(M) the distribution becomes single peaked and the lognormal does also fit an important 
fraction of the low density part. The widths of the lognormals, plotted as a function of 
the Mach number in Figure^ {solid line), are systematically larger than the corresponding 
values found by Burkhart & Lazarian (2012) for the isothermal case {dotted line), but they 
follow a qualitatively similar behavior with M. In fact this behavior can be fit by a function 
of the form (|2]) with A-£ = 0.084 and 6s = 12.5 {dashed red line) obtaining an error in the 
fit of 2.17%. This set of parameters is however not unique. As an example we also display 
the curve for A-^ = 0.081 and 6s = 14.29 {dashed blue line) which has an error of 2.22%. 
Considering the fact that errors in fiting lognormals are greater than the difference between 
the errors resulting from the two previous fits, we can consider them as being equivalent. 
Note that in Figure [6] the Mach number is the rms value at the mean temperature Mrms- We 
choose this value because the average local Mach is dominated by the warm gas (T > 6100K), 
whose volume fraction is between ~ 30%, for very turbulent simulations, and ~ 80%, for 
the smaller value of M, which is very large when compared with the one of the cold gas 
(T < 310 K) which goes from ~ 2 to 8%. 
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4. 



Discussion 



4.1. 



Diffuse Gas 



The approximate lognormal shape of the diffuse gas density PDF is consistent with 
observations reported by Berkhuijsen & Fletcher (2008) for galactic low density neutral gas. 
However, the fact that the lognormal fit fails at very low densities is also consistent with 
predictions by Passot & Vazquez- Semadeni (1998) for polytropic flows with 7 < 1, according 
to which the PDF is expected to decrease faster than a lognormal at low densities. 

The behavior of as with M, that grows faster than in the isothermal case, shows the 
effects of the presence of TI. In fact as discussed in Vazquez- Semadeni et al. (2003), the 
increase of M has two effects on the development of TI. First, as the turbulent crossing time 
decreases it becomes shorter than the growth time for linear TI, which is the time for the 
gas to form condensations, increasing the fraction of unstable gas (e.g. Gazol et al 2001) 
and producing a larger drift from the thermal equilibrium. On the other hand, the presence 
of velocity fluctuations generate adiabatic density perturbations which are linearly unstable 
only when the cooling time is shorter than the dynamical time (which for the supersonic 
case is the turbulent crossing time), i.e. only at large scales . However, in the nonlinear 
regime, reached as M increases, due to the local increase of density in forcing generated 
compressions, the cooling time can locally decrease implying that even initially small scale 
fluctuations can become thermally unstable. The consequence of this nonlinear development 
of TI is an enhancement of density contrast that could lead the width of the density PDF 
to grow faster with M than in the isothermal case. 



The results we obtained for the high density gas reconcile the theoretical prediction from 
Passot & Vazquez- Semadeni (1998) with previous numerical results for thermally bistable 
flows (Audit & Hennebelle 2010). Following the former, for a cooling function adapted 
to describe the dense atomic gas, having an effective polytropic index 7 < 1 in thermal 
equilibrium, the density PDF is expected to develop a power-law at high values. On the 
other hand Audit & Hennebelle (2010) report that the dense gas density PDF resulting from 
thermally bistable colliding flows seems to be better described by a lognormal than by a 
power-law. From results presented in §3. H it is clear that the dense gas PDF can behave as a 
power-law or as a lognormal depending on the M value and more specifically on the amount 
of dense gas out of thermal equilibrium. The transition between these behaviors could be 
due to the fact that the higher is the amount of gas out of thermal equilibrium, the lesser is 



4.2. 



The dense Gas 
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adequate a polytropic equation of state as a thermodynamic description of the gas. In fact, 
Passot & Vazquez-Semadeni (1998) suggest that the development of a power-law for 7 7^ 1 
is the consequence of density jumps depending the local density. In our simulations the 
dense gas drifting away from thermal equilibrium is in fact the consequence of the density 
fluctuations being determined by velocity fluctuations. A remarkable difference between 
the power-laws we obtain at high densities for low M simulations and the predictions from 
Passot & Vazquez-Semadeni (1998) is the logarithmic slope behavior. We get slopes of —4.43 
and —2.41 for (M) = 0.28 and (M) = 0.73, respectively. This implies that the power-laws 
resulting from our simulations are stepper than the power-law they predict for 7 = 0.5 which 
is expected to have a logarithmic slope of —1.2. Note however that this value corresponds 
to the large M limit. 

The development of a lognormal at high densities in high M simulations suggests that 
in some cases the use of this kind of distribution as initial condition in the molecular cloud 
formation process represents an adequate choice. Nevertheless, for our simulations the re- 
lationship (Jg-M is very different from equation ([3]), implying that in early times during the 
cloud formation process the density distribution can not be related with the dynamical state 
of the gas through the isothermal version of this equation. In particular, for a given Mach 
number we find a narrower distribution. It is important to note that the presence of self 
gravity, which is not included in our models, could have noticeable effects on the distribution 
of dense gas even in these early phase. This problem is going to be addressed in a future 
work. 



4.3. The column density 

The dense branch of the column density PDF is well described by a lognormal regardless 
of the M value. This is partially consistent with recent observational and numerical works 
suggesting that in clouds where the effects of gravity are not dominant in determining the 
cloud structure, the column density PDF has a lognormal shape (Kainulainen et al. 2011, 
Ballesteros-Paredes et al. 2011). In those works however, turbulence is invoqued as the main 
agent in shaping the column density PDF. From our results it is possible to suggest that even 
when low levels of turbulence are present, the cooling properties of the gas could produce 
a lognormal column density PDF. In fact, Heitsch et al. (2008) investigate the expected 
timescales of the dynamical and thermal instabilities leading to the rapid fragmentation 
of gas swept up by large-scale flows and compare them with global gravitational collapse 
timescales. They identify parameter regimes in gas density, temperature, and spatial scale 
within which a given instability dominate, finding that the thermally dominated parameter 
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regime has a remarkable large extent due to the fact that outside the strictly thermally 
unstable regions, cooling could still being the dominating agent leading to fragmentation in 
the presence of an external (in this case ram) pressure. 

On the other hand we find that the width of the column density distribution is much 
larger than in the isothermal case studied by other authors for MHD flows (Kowal, Lazarian 
& Brcsnyak 2007; Burkhart & Lazarian 2012). This is a natural consequence of the large 
density dynamical range produced by the presence of thermal instability. The resulting 
(^inCs/Sigmao)-^ relationship has to include the fact that the density contrast is large even 
for low Mach numbers requiring a very rapid grow that implies parameter values completely 
different from the ones reported in the isothermal case, namely As = 0.11 and 6s = 1/3. 

4.4. Numerical Issues 

Several numerical factors can affect in some way the results obtained in the present 
work. 

4.4- i- Resolution 

We have performed some higher resolution simulations in order to see the effects on the 
density PDF. For similar Mach numbers we find that the main consequence of increasing the 
resolution to 1024^ seems a deviation of the high density tail in moderate Mach simulations. 
In particular, for (M) ~ 1.9, the distribution falls faster than a lognormal for high n. 
Unfortunately we do not have enough snapshots to quantify this effect and to measure a in 
high resolution simulations. Although our PDFs could be not fully converged and higher 
resolution simulations could potentially lead to results quantitatively different from those 
presented in previous sections, the facts that the diffuse gas PDF can be well described by 
a lognormal and that the dense gas PDF is well described by this kind of distribution only 
for large enough values of M, does not seem to change with resolution. 

4.4.2. Model 

Even if our model allows both, the study of diffuse gas distribution and a direct com- 
parison with isothermal numerical experiments reported by other authors, there are some 
choices which can potentially affect our results. The periodic boundary conditions maintain 
a fixed amount of gas in the box and this fact can artificially regulate the gas segregation. 
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Also the particular choice of the cooling function can affect the gas segregation. This func- 
tion could change because of physical reasons such as abundance variations, heating rate 
variations, or additional cooling process. Finally, additional physics such as the presence of 
self gravity and magnetic fields could also modify the quantitative behavior of the density 
distribution. 



4.4- 3- Forcing 

Federrath et al. (2008) and Federrath et al. (2010) found that for isothermal gas the 
width of the PDF depends not only on the rms Mach number but also on the relative degree 
of compressible and solenoidal modes in the turbulence forcing, with 6 = 1/3 appropriate 
for purely solenoidal and b = 1 for purely compressive forcing. For the thermally bistable 
case the effects of forcing could be more complex due to the interplay between TI and 
turbulence. Although simulations with different kind of forcing have been presented by 
Seifried at al. (2011), the energy transfer between the solenoidal and the compressive modes 
for turbulent bistable flows have not been addressed. As shown in Gazol & Kim (2010) for 
purely solenoidal forcing, the presence of TI does significantly affect the density as well as 
the velocity power spectrum. In fact, it is well known that the development of TI produce 
turbulent motions with typical velocities of the order of tenths of km s^^ (Kritsuk & Norman 
2002; Piontck & Ostriker 2004) but, as stated earlier, the presence of turbulence does in turn 
modify the development of TI. A detailed analysis of the effects of forcing comparing also 
the resulting velocity power spectrum and the density weighted velocity power spectrum for 
the compressible as well as for the solenoidal modes, should be done in order to address this 
problem, but it is out of the scope of the present work. 

5. Summary and Conclusions 

In the present work we have presented numerical experiments showing that: 

1. At low densities the volume density PDF resulting from turbulent thermally bistable 
flows can be well described by a lognormal distribution whose width increases with the 
Mach number. 

2. A relationship between the width of the distribution and the average local Mach number 
can be found, however this relationship is not the same as in the isothermal case. 
The value of the parameter b, included in the isothermal case, is for our simulations 
surprisingly close to the one obtained for purely solenoidal forcing in isothermal gas 
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but the non linear development of TI, producing a more efficient rarefaction of gas, 
causes a faster growth of the distribution width as M is increased. The consequence of 
this enhanced growth in the mathematical form of the ag — M relation is the presence 
of a scale factor distinct from 1. 

3. At high densities the form of the volume density PDF depends on the value of the 
Mach number. For simulations with transonic or weakly supersonic average velocities 
in dense gas the distribution is a power-law, while in the presence of highly supersonic 
velocities the distribution becomes lognormal. 

4. For simulations which develop a high density distribution with a lognormal shape, the 
witdth of the distribution is smaller than in the isothermal case and grows slower with 
the Mach number. 

5. At high densities the column density PDF resulting from our simulations can be de- 
scribed by a lognormal for all the Mach numbers we consider. As M increases the 
density range where the lognormal fit is adapted expands and the lognormal becomes 
wider. 

6. The width of the column density distribution resulting from our simulations is system- 
atically larger than the width obtained in the isothermal case. A relationship between 
the width of the column density distribution and the rms Mach number at the mean 
temperature can be found. This relationship has the same form as the one reported 
in the literature for the isothermal case, but the parameters resulting from our fit are 
very different. 

From these results it is clear that when studying the diffuse and/or the dense atomic 
interstellar gas in order to relate the density structure, and in particular the width of its 
PDFs, with the dynamical state of the gas, characterized by the Mach number, the use of 
results obtained from isothermal turbulent flows is not an adequate choice. Specific relations 
between as (o"i„(i;/So)) ^ for thermally bistable fiows should be taken into account in 
any observational or theoretical work using the density (column density) PDF as a measure 
of the Mach number. The relationships obtained in the present work could be affected by 
the inclusion of additional physics such as self gravity, magnetic fields or variations on the 
cooling function due to variations of the heating rate and the gas abundances. 

The work of A. G. was partially supported by UNAM-DGAPA grant IN106511 Some 
of the numerical simulations were performed at the cluster Platform 4000 (KanBalam) at 
DGSCA, UNAM. 
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Fig. 1. — Density PDF for different values of (M). Dashed and dotted lines represent 
lognormal fits to the high and low density parts of the distribution, respectively; while 
dashed-dotted lines in low (M) PDFs correspond to power-law fits with logarithmic slopes 
of -4.43 and -2.41 for (M) = 0.28 and (M) = 0.73, respectively. 
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Fig. 2. — Lognormal widths as a function of for fits corresponding to low density gas 
in PDFs displayed in Figure [H Dashed black lines plot (jg = (ln(l + b'^M'^)y^'^ for b = 1/2 
and 1/3, whereas dashed red line corresponds to the fit described in the text. Mach number 
values are average values of the local Mach number for the diffuse gas. 
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Fig. 3. — Lognormal widths as a function of Mc for fits corresponding to the high density 
gas and (M) > 1.86 in PDFs displayed in Figured] 
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Fig. 4. — Effective polytropic index 7 as a function of the corresponding averaged local Mach 
number for the whole simulation {solid line, (M)), the diffuse gas {dashed line, M^) and the 
dense gas {dotted line, M^). 



Fig. 5. — Thermal pressure-density relation for simulations with (M) = 0.28 (left) and 1.86 
(right). In each panel the black curve corresponds to the thermal equilibrium implied by 
the cooling function and contours are placed at 10% (violet), 30% (blue), 50% (green), 70% 
(orange), and 90% (red) of the logarithm of the maximum value of the two dimensional 
histogram. 
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Fig. 6. — Column density PDF for different values of (M). Dotted lines represent lognormal 
fits. 
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Fig. 7. — Lognormal widths as a function of Mrms for fits corresponding to column density 
distributions displayed in Figure O In this plot Mach number values are rms values at the 
mean temperature. Dashed black /mes plot equation ([2]) for = 0.11 and 6s = 1/3, whereas 
dashed red and dashed blue lines correspond to the fits described in the text. Mach number 
values are the rms values at the mean temperature Mrms- 



